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ABSTRACT 

We re-examine claims for redshift evolution in black hole-bulge scaling relations based on lensed quasars. 
In particular, we refine the black hole mass estimates using measurements of Balmer lines from near-infrared 
spectroscopy obtained with Triplespec at Apache Point Observatory. In support of previous work, we find a large 
scatter between Balmer and UV line widths, both Mg II AA2796,2803 and C IV AA1548, 1550. There is tentative 
evidence that C III] A 1909, despite being a blend of multiple transitions, may correlate well with Mg U, although 
a larger sample is needed for a real calibration. Most importantly, we find no systematic changes in the estimated 
BH masses for the lensed sample based on Balmer lines, providing additional support to the interpretation that 
black holes were overly massive compared to their host galaxies at high redshift. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: Seyfert 



1 . EVOLUTION IN BLACK HOLE-BULGE RELATIONS 

Locally, we observe tight correlations between the properties 
of bulge-dominated galaxies and the masses of their central 
supermassive black ho l es (BHs; e.g., Kormendv & Richstond 
[19951: iTremaine et all 12001 iGultekin et al] l2009h . The 
mechanisms that establish and maintain these relations 
are uncertain, despite innumerable suggestions i n the 
literature (e.g.. ISilk & Red 119981: iMurravetalJ [20051 



Miralda-Escude & Kollmeier 12005F Hop kins et alj|2006t iPen 



20071) . [n principle, evaluating the demographics of nuclear 
BHs as a function of redshift should observationally constrain 
the processes that lead to the tight scaling relations observed 
today. Unfortunately, it is currently prohibitive to obtain dy- 
namical BH masses for systems beyond tens of Mpc. Thus, 
estimates of BH mass at large distance are necessarily based 
on very in direct methods link ed to accretion processes in active 
BHs (e. g.. I Vestergaardll2002h . 

Several studies have used active galaxies to probe evolution 
in B H-bulge scaling relations, probing redshifts from .4 < z < 
0.6 dWoo et al]2006tlTreu et alj2007tlWoo et alj2008l) and 1< 
z < 4 dPeng et alJ|2006allUISalviander et alJl2007tlShields etal] 
2006: lHdl2007tlJahnke et ^l.ll2009t iMcLeod & Bechtoldl l2009) 
all the way to z = 6.4 (IWalter et al.l2004l) . Generally speaking, a 
wide variety of observations suggest that BH-bulge relations do 
evolve with redshift, in the sense that the r atio of BH to bulge 
mass was higher at ear ly times (but see also Shi elds et al.| [2003; 
lAlexanderetai] 120051) . although we have observational con- 
straints only for the most massive systems (Mbh> 10 8 Mq) at 
high redshift. This counterintuitive result has stimulated vigor- 
ous discussion both abou t the ramifications for the coevolution 
of BHs and bulges (e.g., iRobertson etal] 120061: ICrotonll2006l) 
and about whether there are b uilt-in biases in the measurement 
techniques dLauer et al. 2007). 

Unfortunately, even apart from potential population biases, 
our interpretation of the observations are prone to significant 
uncertainty. On the one hand, it is only possible to ob- 



tain BH mass estimates at cosmological distances using ac- 
tive galaxies (typically luminous quasars at high redshift). 
Immediately it becomes very challenging to characterize the 
host galaxy properties, when the quasar outshines the under- 
lying galaxy starlight by factors of ~ 10-30 (e.g.. |Peng et al.1 
2006at iKim et alj|2008af) . Apart from direct imaging, some 
groups have used gas me asurements (predominantly CO) to 
obtain dynamical masses dShields et all 120061: iHo et al.1 120081; 
IWalter et all 12004 iRiechers et al.1 120081 120091) . which may o r 
may not provide a reliable tracer of the galaxy mass 

dHc|2007). 

The width of narrow emission lines, particularly [O III] A5007, 
have also been substituted for the galaxy velocity disper- 
sion (e.g.. IShields et al.l l2003l; lBorosonll2003l: ISalviander etail 
2007; Gaskell 2009). While there is a strong correlation be- 
tween stellar a nd gaseous velocity dispersion in low-lum i nosity 
sources (e.g iHeckman et al.1 119811 iNelson & Whittle! fl996: 
iGreene & Holl2005al: iHol 120091 there is good r eason to sus- 
pect t ha t it does no t hold at high luminosity (e.g. jGreene et al.1 
2009!). iPeng et all d2006bl P06 hereafter) mitigated the host 
galaxy contrast problem by focusing on lensed quasars. The 
quasars are lensed differently from the underlying (resolved) 
host galaxies, reducing the contrast problem discussed above. 
For that reason, we focus on the lensed quasar sample in this 
paper. 

Daunting as measuring high-redshift galaxy mass and stel- 
lar velocity dispersion may be, particularly in the presence of 
a luminous quasar, the BH mass measurements are equally 
problematic. The techniques are indirect and model-dependent. 
Briefly, active galaxies contain dense gas orbiting at distances 
of light days to months from the central BH that gives rise 
to broad emission lines with widths of thousands of km s . 
By combining a size scale with the line width of the emit- 
ting region, the broad-line region (BLR) gas can be u sed as 
a dynamical tracer of the BH mass re.g.. [Pibail [l980h . Di- 
rect size estimates are obtained by measuring the time lag be- 
tween variability in the continuum and line emission (reverbera- 
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tion mapping : lBlandford & M cKee 1982; Peterson et al]|2004l) . 

Table 1. Sample and Observations 



Galaxy 


RA 


Dec. 


m H 


Obs. Date 


'exp 


SIN„ a 


S/N„ a 


<i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


PMNJO 134-0931 


01:34:35.7 


-09:31:02 


15.3 


2009-07-14 


3240 


295.8 


13.7 


Q0 142- 100 


01:45:16.5 


-09:45:17 


17.6 


2008-10-18 


2880 


63.3 




SDSS0246-0825 


02:46:34.1 


-08:25:36 


18.1 


2008-11-11 


4320 


1 1 1.2 


14.3 


MG0414+0534 


04:14:37.7 


+05:34:44 


15.9 


2008-10-18 


3600 


140. 1 


15.5 


HS0S10+2554 


08:13:31.3 


+25:45:03 


14.8 


2008-10-18 


1800 


393.7 


54.3 


SBS0909+523 


09:13:01.1 


+52:59:28 


14.7 


2008-10-18 


2 1 60 


374.6 


163.7 


FBQ095 1+2635 


09:51:22.6 


+26:35:14 


17.0 


2008-11-11 


2160 


45.0 


19.2 


Q0957+561 


10:01:20.8 


+55:53:49 


15.6 


2008-11-17 


2880 


460.1 


111.1 


J1004+1229 


10:04:24.9 


+12:29:22 


17.8 


2008-11-17 


5400 


125.9 


17.2 


HE1 104-1805 


11:06:33.5 


-18:21:24 


15.9 


2009-03-14 


2880 


1 82.9 


10.1 


PG1 115+080 


11:18:17.0 


+07:45:57 


15.7 


2008-03-24 


1440 


152.8 


31.8 


B 1152+200 


11:55:18.3 


+19:39:42 


15.6 


2009-06-05 


3600 


50.6 


9.2 


H1413+117 


14:15:46.4 


+11:29:41 


15.8 


2009-05-05 


3600 


137.2 


18.8 


B1422+231 


14:24:38.1 


+22:56:00 


14.4 


2009-05-05 


1440 




39.9 


FBQ1633+3134 


16:33:49.0 


+31:34:11 


15.8 


2009-05-11 


1440 


104.0 


13.2 


Q2237+030 


22:40:30.3 


+03:21:28 


15.0 


2009-07-14 


3240 


128.8 


27.9 



Note. — Col. (1): Name. Col. (2): Right Ascension (hrs; J2000). Col. (3): Declination (deg; J2000). 
Col. (4): Observed //-band magnitude (mag). Col. (5): Date of observation. Col. (6): Total on-source 
exposure time (sec). Col. (7): Signal-to-noise ratio summed over the Ha emission line, excluding any 
regions that were excised from the fit. Col. (8): Same as Col. (8) for the H/3 line. 



A tight, empirically determined correlation between BLR radii 
and the l uminosity of the active galaxy (the radius-lu minosity 
relation; iKaspi et al.l 12005b IBentz et alJl2006l l2009al) can be 
used to obtain a pproximate BLR r adii for the quasar population 
in general (e.g.JVestergaard 2002). The current data support a 
slope of Rblr oc L 5 , as is expected if neither the density struc- 
ture of the BLR nor the s hape of the ionizin g continuum de- 
pends on luminosity (e.g.. IBentz et al.l [20061) . Unfortunately, 
the radius-luminosity relation has not been calibrated exten- 
sively for lumino sities greater than Lsioo~ 10 46 erg s -1 (see 
IKaspi et al.ll2007h . 

Because the kinematic structure and inclination of the BLR 
are unknown, the derived virial "masses" (Mbh<x- v 2 Rblr/G) 
have no physically motivated normalization. Recent reverbera- 
tion mapping experiments reveal signatures of inflow, rotation, 
and outflow in indi vidual objects, but such two-dimensional 
maps remain scarce dWelsh & Hornelll99ll ; IBentz et aU l2009c ; 
iDennev et alj l2009blh Our current practice is to compare 
the virial masses with independent estimates of BH mass, 
typically using th e Mbh ~ o~* relation, to derive an average 
scale factor (e.g., Gebhardt et alJ 12000 ; Ferrarese et al. 2001; 



Onkenetal. 2004; Nelson et al. 2004; Greene & Ho 2006; 
She n^tajMgi fact that any^elation is seen betw^n 
virial masses and the mass inferred from the Mbh - c* relation 
is encouraging, but leaves much room for large systematic un- 
certainties (e.g jKrohkfeOQlUCollin et al.ll2006l) . 

At higher redshift, only rest-frame UV spectra are read- 
ily available for large samples of quasars. The scaling rela- 
tions for the UV lines (e.g., Mg II A2800 A; C IV A1550 A) 
include additional layers of uncertainty. Virtually no rever- 
beration mapping has been performed with Mg II, and the 
scali ng relations for this line a re simply scaled to match H/3 
(e.g.. iMcLure & Dunlopl 120041; lOnken & Kollmeierll2008l) . In 
the case of C IV, re verberation mapping has been done (e.g., 
iPeterson et al. 1 12005b . but there are strong reasons to suspect 
that the C IV line width is not dominated by virial motions 
( e.gjGaskelJll982HBaldwin et al.ll 19961; [Richards et al. 2002a; 

eighly & Moord 12004 iBaskin & Laorl 120051; ISulentic et al.l 
2007j IShen et all 12 008b) al though debate continues on this 



point (e.g., [Ve stergaar d& Peterson! 120061; iKellv & Bechtoldl 
12001; Gavigna udetalJl2008l) . For these reasons, virial masses 



based on Balmer lines (preferably Ha; iGreene & Holl20 05b) 
have the most credibility, since these have been directly com- 
pared with alternate estimates of Mbh- We focus specifically 
on obtaining Balmer-based virial masses for the high-redshift 
lensed quasar sample from P06. Our primary goal is to deter- 
mine whether the masses presented in P06 are systematically 
biased by the use of UV line transitions. We follow P06 and as- 
sume a standard cosmology with Ho = 100 h = 70 km s" 1 Mpc" 1 , 
fl m = 0.30, and fl A = 0.70. 

2. OBSERVATIONS AND DATA REDUCTION 

The data presented here were obtained over the course of a 
year using t he newly commis sioned near-infrared spectrograph 
Triplespec ( Wilso net al J I2004T) at Apache Point Observatory 
(Table 1). All objects were observed with a 1.1 x 43" slit and 
Fowler sampling of 8 (the number of non-destructive readouts 
at the begi nning and end of each e xposure designed to minimize 
readnoise; lFowler & Gatlevl 19901) . Triplespec covers a nominal 
wavelength range of 0.95-2.46 /im with R = 3500. Observing 
conditions ranged from clear to partly cloudy, with typical see- 
ing of 9 » l."5. In most cases the slit was positioned at the 
parallactic angle in the middle of the observation, although in a 
couple of cases we positioned the slit to place two quasar im- 
ages in the slit at once. The object was dithered along the slit 
every 180 sec to improve sky subtraction. For each quasar we 
observed a nearby A0V star (10 < H < 6 mag) to serve as flux 
and telluric standard. 

The data were reduced using custom software th at is a modi- 
fied ve rsio n of Spextool and is de scribed in detail in lVacca et al.l 
d2003l) and lCushing et alj d2004l) . Using dome-flat and arc-line 
exposures, the code creates and applies a flat-field and wave- 
length solution. Bias and dark subtraction is accomplished 
through pair-wise differencing of images taken at two slit po- 
sitions, which also removes air glow emission from the atmo- 
sphere, at least to zeroth order. Nonlinearity corrections are ap- 
plied a nd then each pair of spectra are traced and optimally ex- 
tracted dHornel 1986b . including background subtraction. Wave- 
length calibration is applied, and all the spectra of a given 
source are median-combined. Flux calibration is accomplished 
using an A0V star observed at similar time and airmass. This 
same star is used to create a model of the telluric absorption by 
assuming that the A star has an intrinsic spectrum identical to 
that of Vega (see Vacca et al. for details). Prior to correction, 
small wavelength shifts between the A star and the program 
object are derived on an order by order basis using a cross- 
correlation technique. Finally the orders are merged with small 
scale factors applied to properly match the edges of each order, 
and cosmetic data clipping is done. The software also generates 
an error array. Heliocentric corrections are calculated using the 
IRAF task bevcorr. 

The resulting S/N for each object across the Ha and H(3 lines 
are shown in Table 1. The relative flux calibration is reason- 
able and the resulting spectra have smooth, power-law continua. 
However, given the variable clouds that plagued many of our 
observations, we suspect that the overall flux calibration scale 
is not reliable. For instance, the flux scales for two observations 
of B 1152+200 differ by a factor of three. Our results are not im- 
pacted by these problems, however, since we measure intrinsic 
luminosities from broad-band photometry combined with a lens 
model (P06 and see below). 

2.1. Rest-frame Ultraviolet Spectra from SDSS 
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FIG. 1. — Example fits to the continuum (left), H/3 (middle), and Hct (right) lines from the Triplespec data. We show the data (solid histogram), the total model 
(thin solid), the broad- (dashed) and narrow-line (dotted) model components, and residuals below (thin solid histogram). Data are plotted with an arbitrary scale in 
f\. The rest of the sample are plotted in the Appendix. 
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FIG. 2. — Example fits to the continuum (left), C iv or C m] (middle), and Mg n (right) lines from the SDSS data, shown in units of 10~ 15 erg s -1 cm -2 A -1 . 
We show the data (solid histogram), the total model (thin solid), the broad- (dashed) and narrow-line (dotted) model components, and residuals below (thin solid 
histogram). The dashed vertical lines denote a masked region of the spectrum that is heavily absorbed in the original spectrum. Note that in the C m] profile 
(top-middle panel) we show only the broad component of the C in] line itself. 
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In addition to the rest-frame optical spectra obtained with 
Triplespec, we also utilize observed optical (rest-frame ultravi- 
olet; UV) spectra from the the Sloan Di gital Sky Survey (SDSS; 
lYork et al]|2000l ; lAbazaiian et al.| [2009). These t argets were se- 
lected from the SDSS photometry as quasars (Ric hards et al.l 
2002b). Galaxies with SDSS observations have tabulated UV 
slopes in Table 2. 

3. CONTINUUM AND LINE MEASUREMENTS 



3.1. The Fits 

Quasar continua are generally well-fit by a power law. How- 
ever, superimposed on this smooth continuum is a "pseudocon- 
tinuum" of broad Fe II multiplet emission that effectively litters 
the entire optical/UV region of the spectrum. Although much 
progress has been made deriving theoretical Fe II spectra (e.g., 
IVerner et al. 2004), the prospect of fitting > 800 transitions in- 
dividually is a daunting task. Rather, it is common practice 
to derive an Fe II template from a high S/N observation of an 
active galaxy whose broad emission lines are intrinsically nar- 
row (^5 1000 k m s" 1 ), typically 1 Zw I Here we use the opti- 
cal template of iBoroson & Greed (1 19921) and the UV template 
pr esented in Salviander et al. (2007) that was derived from that 
of [vistergaard & Wilkesl ifeOOll) . but has a theoretical templat e 
pasted in beneath the Mg II line from lSigut & Pradhanl d2003h . 
Vestergaard & Wilkes also provide a separate Fe III template at 
k 1900 A that we utilize in at tempting to fit theC III] A 1909 A 
line. In addition, we follow Dietrich et I] (|2002) and model 
the "little blue bump' ' (lGrandilll982j[ as a combination of 
higher-order Balmer emission lines dStorev & Hummeilll995h 
and bound-free Balmer continuum emission with a temperature 
T e = 15000 K, electron density n P = 10 8 -10 1 cm" 3 , and optical 
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depth 0.1 < r v < 2 (following lGrandilll982l) . 

In practice, we fit the power-law continuum, Fe II emission, 
and Balmer continuum simultaneously, using emission-line- 
free windows (see Figure 1). The fits include a normalization 
and slope for the (single) power-law component, a width, shift, 
and amplitude for the Fe template, and a density and optical 
depth for the Balmer continuum. In the UV, we find the best 
results fixing the Fe II broadening to that of H/3, while in the 
optical we find negligible difference in allowing the broaden- 
ing to be free. We subtract the continuum model and create a 
line-only spectrum. 

We fit the emission line spectra with multi-component Gaus- 
sians. These components are used only to create a high-fidelity 
noise-free match to the emission lines; we do not ascribe phys- 
ical meaning to them. We first fit the H/3+[0 III] AA4959,5007 
region. The narrow [O III] lines are fit with up to three Gaus- 
sians each, with the relative wavelengths tied to laboratory val- 
ues and the line intensities constrained to have a flux ratio of 
1:3 (T able 2; see details in IGreene" & Holl2005at iGreene et al.1 
2009). The width of the narrow H/3 line (fit with a single Gaus- 
sian) is tied to that of [O III] whenever possible (Table 3fl 
Finally, broad H/3 is fit with up to four (typically two) broad 
components. We also find the need to impose a lower limit 
of 1000 km s -1 on the broad components so that they do not 
erroneously fit a narrow-line component. We note that in these 
bright quasars the flux contribution from the narrow lines is typ- 
ically small. 

'SBS0909+523 is the one exception. The narrow line component is very strong in this object, even in the UV lines (see Figs. 1 & 2) and is generally narrower in 
other transitions than in the [O III] line. Interestingly, the broad lines appear to be redshifted compared to the narrow lines in this object. It would be worth attempting 
integral-field spectroscopy to investigate any spatial offsets corresponding to the observed velocity offset. 



2000 4000 6000 8000 
FWHM Ha (km s" 1 ) 

FIG. 3. — Comparison between the line widths derived from Ha and H/3 
in the Triplespec data. H/3 lines with incomplete spectral coverage are noted 
with a cross. The one-to-one line (solid) is shown for reference only. As 
expected (e.g., Greene & Ho 2005b) H/3 is generally broader than Ha. Ex- 
cluding the compromised H/3 fits, we find (FWHM Ha /FWHM H /3 > = 0.9 ±0.2 
(d ashed line), wh i ch is c onsistent with the 7% difference in line width found 
by |Greene & Ho ( 2005b), albeit with considerable scatter. We note that 
MG0414+0534 has a very large and very unreliable H/3 line that is off the 
scale of this plot. 



Nevertheless, the treatment of the narrow lines is a source of un- 
certainty in the FWHM measurements, particularly in the UV. 
As described below, we therefore perform an additional fit with 
no narrow component and fold the difference into our total error 
budget. 

Our a pproach is similar in the case of the H a+[N II] complex 
(see also lHo et aUllWk IGreene & Holl2004l) . Here we fix the 
narrow Ha and [N II] to the [O III] line width, fix the relative 
wavelengths to laboratory values, and fix the relative strengths 
of the [N II] AA6548,6584 A to 1:3. In cases where there is no 
[O III] line, we fix the narrow-line width to 500 km s" 1 . Gen- 
erally in these cases the narrow components are too weak to be 
fit independently. Uncertainties in this procedure are estimated 
using an alternate fit with no narrow-line components. We fit 
the broad Ha with as many as four Gaussians, but do not tie 
them to the H/3 profile in any way. All of the Triplespec spectra 
and fits are shown in Figure 1 and the Appendix. 

Our fits to the UV lines proceed in a similar fashion, with 
each broad line modeled as the sum of up to four Gaussians. 
We still choose to tie the narrow-line components to the width 
of [O III] as above, although we note that the proper treatment 
of the narrow co mponent of C IV remains a matter of debate 
in the literature dBaskin & La or 2005 : IVestergaard & Peterson! 
120061: iKellv & Bechtoldl 120071: IShen et alj|2008bl) . Again, in 
cases without available [O III] fits we have simply chosen a 
representative width of 500 km s _1 , and again we perform a 
second, narrow-line-free, fit. Finally, we include a linear con- 
tinuum component to remove residual continuum errors (Fig. 
2). We note that the quality of the Mg II fit does depend on 
the Fe II subtraction, and our model for the Fe II continuum is 
poorly constrained directly beneath the Mg II line. 
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Table 2. Spectral Measurements 



Name 


F< Wo[III] 


FWc n 
vv Fell 


vv H/3.n 


^"11/3.1 


FWu 




FWu , 
1 : vv Hct.t 


C^opt 


CKUV 










V+) 


C5~l 




("71 


K°) 


\") 






PMNJ01 34-0931 


12 


45 


2.1 


68 


70. 


9.6 


491 


-0.4 






Q0142-100 




11 




42 


<0.8 


<0.5 


212 


-2.2 


-1.5 


-1.8 


SDSS0246-0825 


6 


58 


1.5 


50 


<0.6 


4.8 


212 


-2.4 






MG0414+0534 




163 


0.7 


123 


3.6 


1.3 


340 


0.10 




1.0 


HS08 10+2554 


10 


76 


<0.5 


86 


<0.7 


<0.5 


377 


-1.3 


-1.8 




SBS0909+523 


63 





7.5 


119 


<1.2 


16. 


526 


-1.3 


-0.9 


-1.0 


FBQ0951+2635 


43 


157 


<1.2 


128 


<2.9 


9.5 


340 


-1.4 


-1.0 


-1.7 


Q0957+561 


14 


25 


1.0 


43 


<0.5 


3.4 


158 


-1.5 


-1.8 


-1.8 


J1004+1229 




60 


9.4 


76 


<1.1 


0.9 


327 


-0.6 




-1.5 


HE1 104-1805 


9 


31 


1.2 


64 


12. 


8.7 


230 


-2.5 




-2.0 


PG1 115+080 




48 




50 


<0.7 


<0.5 


289 


-1.8 


-1.8 


-2.0 


B 1152+200 










37. 


15. 


744 


4.10 


0.1 




H1413+117 


39 


20 


10. 


75 


<1.2 


40. 


314 


-1.7 


-0.3 


-0.9 


B 1422+231 


26 





9.6 


60 








-1.5 




-1.4 


FBQ1633+3134 


29 


28 


6.4 


62 


<0.9 


<0.5 


325 


-1.1 


-2.1 


-1.8 


Q2237+030 


12 


10 


0.7 


30 


11. 


5.2 


114 


-1.5 







Note. — Col. (1): Name. Col. (2): Equivalent width (EW; A) of the [O m] A5007 line. Here and in the following entries 
the EWs are measured from the multi-Gaussian fits. Col. (3): EW (A) of the Fe n line measured between 4434A and 4684A, 
following Boroson & Green (1992). Col. (4): EW (A) of the narrow component of H/3. Col. (5): EW (A) of the total H/3 line. 
Col. (6): EW (A) of the narrow component of Ha. Col. (7): EW (A) of [N n] A6584 A. Col. (8): EW (A) of the total Ha line. 
Col. (9): Slope of the rest-frame optical continuum measured from the Triplespec spectra, fx tx X a °p<. Col. (10): Slope of the 
rest-frame UV continuum measured from the SDSS spectra. Col. (11): Slope measured by P06 from broad-band photometry (see 
text). 



Table 3. Line-width Measurements 



Name 


z 


FWHM [0 ni] 


FWHMhq 


FWHM H/ 3 


FWHMm.ii 


FWHMcni] 


FWHMciv 


Ai5100 


log M BH ,P06 


log MBH.Ha 


logM BH ,H/3 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(ID 


(12) 


PMNJ0 134-0931 


2.22 


1410 


5.1±0.5 


4.2± 1.4 
















Q0142-100 


2.7 




3.8± 0.3 


2.7± 0.6 




6.0± 0.2 


6.5± 1.3 


46.42 


9.4 


9.4± 0.3 


9.1 ±0.4 


SDSS0246-0825 


1.69 


540 


2.5± 0.2 


2.5± 0.6 


2.9±0.1 


4.5± 2.8 


5.0± 0.2 










MG0414+0534 


2.6 




5.3± 0.3 


10.8±2.3 








45.73 


9.3 


9.4± 0.3 




HS0810+2554 


1.51 


490 


3.8± 0.0 


4.4± 1.1 


3.6±0.1 


4.6± 0.4 


3.9± 0.1 










SBS0909+523 


1.38 


410 


3.1±0.6 


4.1± 1.6 


3.5± 0.5 


4.1± 1.1 




45.96 


9.6 


9.0± 0.3 


9.2± 0.5 


FBQ0951+2635 


1.25 


290 


6.3± 0.4 


7.6± 0.8 


6.0± 0.3 


7.2± 0.8 




45.48 


8.9 


9.4± 0.3 


9.5±0.3 


Q0957+561 


1.41 


570 


3.0± 0.2 


3.3± 0.9 


3.3± 0.2 


5.2± 0.5 




46.25 


9.3 


9.1 ±0.3 


9.2± 0.4 


J1004+1229 


2.7 




3.8±0.1 


5.9± 1.5 








46.12 


9.3 


9.3± 0.3 


9.6± 0.4 


HE1 104-1805 


2.32 


970 


4.7± 0.2 


3.8± 0.9 








45.59 


9.4 


9.2± 0.3 


8.9± 0.4 


PG1 115+080 


1.7 




4.0±0.1 


4.4± 0.2 


4.0± 0.2 


5.7± 3.1 


6.7± 1.0 


45.15 


9.0 


8.8± 0.2 


8.8± 0.3 


B 1152+200 


1.0 




7.1± 2.6 








3.3± 0.6 










H1413+117 


2.55 


1650 


5.3± 0.8 


6.7± 1.9 




8.0± 3.2 


3.3± 0.6 


45.67 


8.4 


9.3± 0.3 


9.5±0.3 


B 1422+231 


3.63 


2720 




6.1±2.2 








46.65 


9.7 




9.9± 0.4 


FBQ 1633+3 134 


1.52 


1720 


4.1 ±0.7 


4.6± 0.9 


4.1 ±0.5 


6.4± 0.6 


3.8± 0.7 


45.81 


9.2 


9.2± 0.3 


9.2± 0.6 


Q2237+030 


1.70 


1070 


4.8± 0.6 


3.8± 1.4 

















Note. — Col. (1): Name. Col. (2): Redshift. Cases with three significant digits were measured from [O in] while cases with two come from the Castles website 
(http://www.cfa.harvard.edu/castles/). See P06 for references. Col. (3): FWHM of [O in] A5007 A line (km s _1 ), with the instrumental resolution of cr* £=: 36 km s -1 sub- 
tracted in quadrature. Col. (4): FWHMn a (10 3 km s" 1 ) measured from the Triplespec data. Col. (5): FWHM H| g (10 3 km s" 1 ) measured from the Triplespec data. Asterisks denote 
objects with incomplete spectral coverage of the H/3 line. Col. (6): FWHMMgii (10 3 km s" 1 ) measured from the SDSS spectrum. Col. (7): FWHMcni] (10 3 km s" 1 ) measured 
from the SDSS spectrum. Col. (8): FWHMciv (10 3 km s -1 ) measured from the SDSS spectrum. Col. (9): Monochromatic luminosity (erg s" 1 ) at A = 5100 A as measured from 
photometry (P06). Col. (10): Log (Mqu/Mq) from P06; typically based on the C iv line. We do not quote errors in the BH masses, since systematic uncertainties dominate. We 
typically quote a factor of four error as derived by Vestergaard & Peterson 2006. Col. (11): Log (Mbh/^o) calculated using Equation 1 and the FWHMH tt from Col. (4) combined 
with the luminosity from Col. (9). Note that the errors in the Ha masses shown here include only measurement uncertainties. Col. (12): Log (M^\\IM(-, ) based on H/3. In this case 
Equation 2 and FWHM H| g from Col. (5) are combined with the luminosity from Col. (9). 




4000 6000 8000 
FWHM Mg „ (km s" 1 ) 

FIG. 4. — Comparison between the line widths derived from Mg 
ii and C m] (black open circles) or C iv (blue open triangles) from the 
SDSS data. The one-to-one line (solid) is shown for reference only. 
We find {FWHMMeil/FWHMcnij) = 0.7 ±0.1 (black dashed line) and 
(FWHMMgii/FWHMciv) = 0.8 ± 0.2 (blue clot-dashed line). 
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FIG. 5. — Comparison between the line widths derived from Ha and Mg 
in or a scaled C m] using a ratio of 0.7 (Fig. 4) when necessary (red open 
squares). As above, we show the one-to-one relation as a solid line. Within 
the scatter, the two measures agree, with (FWHMHa/FWHMMgii) = 1.0±0.2. 



However, numerous authors have now shown that the uncer- 
tainty in F WHM incurred by differing Fe II mode ls is small 
(< 0.02 dex; ISalviander et alJl2007HFine et al.ll2008l) . 

In addition to the standard Mg II and C IV lines, we 
have attempted to model the C III] A 1909 transition, since it 
seems worth investigating every possible transition in the UV 
[and at least in one object, reverberation mapping with this 
line yielded a ma s s that is consistent with other transitions; 
Onken & Petersonl d2002f) 1 . There are good reasons to avoid 
C III], including blen ding with S III] A 189 2, Al III A1857, and 
Fe III multiplets (e.g.. lDietrich et alj|2003) . To minimize these 
degeneracies, each of the three lines (C III], S III], and Al III) is 
modeled as a single Gaussian with the same width. The relative 
centroids are fixed to laboratory values and the strengths are left 
free. We include a narrow component only for C III], and it is 
fixed to the [O III] width as abo ve. Although we also attempted 
to include the Fe III multiplets dVestergaard & W ilkes 200JJ) in 
the fit, keeping the width constrained to that of the Fe II mul- 
tiplets, the results were not well-constrained. It can be seen in 
Table 3 that, unexpectedly, in three out of six cases FWHMcm] 
is larger than FWHMqv- However, these are systems with nar- 
row or broad absorption-line systems, which make the C IV fits 
particularly uncertain. 

After fitting, our velocity measure of choice is a non- 
parametric full-width at half maximum (FWHM) derived from 
our multi-Gaussian fits. We are aware that many authors advo- 
cate the use of the line dispersion rather than the FWHM (e.g ., 
iPeterson et a"i1l2004t lOnken et al.ll200l: iDennev et alj|2009ah . 
We do not believe there to be a strong argument in favor of 
one or the other line-width me asurement at the current time (al- 
though see ICollin et"aI1 2006), but we do know that our tech- 
nique is robust in the presence of noisy spectra, even when fit- 
ting thousan ds of spectra. We dire ct the interested reader to the 
Appendix of iGreene & Hoi d2007al) for a detailed explication of 
our reasoning, but see also lDennev et alj d2009al) . 



3.2. Linewidth Comparisons 

In this subsection we present figures comparing various esti- 
mates of line width derived from different elemental transitions 
(Figs. 3-6). It is clear from these figures that measurement un- 
certainties alone lead to a significant amount of scatter and that 
the total dynamic rang e in linewidth fo r any transition is o nly a 
factor of < three (e.g jFine et al] |2008: She n et al.ll2008bl) . 

For each comparison figure, we have calculated the mean and 
standard deviation in the ratio of the two lines (not the error in 
the mean). These numbers are quoted in the figure captions. 
We have also calculated non-parametric correlation coefficients 
for each line pair. However, presumably because the samples 
are so small, none of the correlations are formally significant. 

We cannot come to any strong conclusions based on these 
comparisons, due to both the small sample size and the large 
measurement errors. We confirm that H/3 i s typically broader 
than Ha (Fig. 3; e.g.. IGreene & Ho1l2005bl , and references 
therein). It has long been known that high-ionization lines 
such as C IV are typically broadened and blu e-shifted com- 



pared to, e.g., Mg II a nd the Balmer lines (e.g., GaskeH| [l982; 
lOsterbrock & Shuded[l982t iBaldwin et ai]|l996l) . While our 
observations are consistent with that trend (Fig. 4), we can- 
not say much more than that. Over this limited range, we do 
not see strong evidence of a correlation between the widths 
of the Balmer lines and C IV, but to a large degree the scat- 
ter is driven by the difficulty in measuring a reliable width for 
the broad-absorption system in H1413+117. Furthermore, the 
(few) C III] measurements we have seem to be correlated with 
the Mg II line widths (Fig. 4). We believe C III] merits inves- 
tigation with a much larger sample. Finally, we tend to derive 
somewhat broader line widths than P06 from the UV lines, pre- 
sumably because we have subtracted the Fe II emission, which 
P06 could not do (Fig. 6). In a couple of cases the difference in 
narrow-line treatment also plays a role. We plan to revisit the 
UV-based mass estimates for the entire sample using modern 
spectroscopy (e.g., from SDSS) in future work. 
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FIG. 6. — Comparison between our measured UV line widths from Mg 
ii (black open circles) or C iv (blue open triangles) and those published in 
P06. Again, the scatter is quite large (FWHM Mg ii/FWHMpo6> = 1.5 ±0.5, 
and again the unity relation is shown as a solid line. Our tendency to derive 
broader line widths than P06 derives from our continuum subtraction and treat- 
ment of the narrow lines. H1413+117 is highlighted as a filled square because 
the broad absorption feature in the C iv line makes fitting difficult. 

3.3. Mass Estimates 

The primary goal of our study is to remove any potential 
systematic bias in the C IV-based masses presented in P06 by 
calculating Balmer-based virial masses. In calculating Ha- 
based BH mass es, we start with the radius-luminosity relation 
of lBentzetaTl(l2009al) . Since this relation was calibrated us- 
ing H/3, we then convert FWHMh0 to FWHMh q using the re- 
lation derived in iGreene &Hol d2O05b), 

which both increases 
the scaling with FWHMh q to the 2.06 power and slightly 
changes the prefactor. Finally, for consistency with P06, we 
assume a scaling factor that is 1.8 times higher than the as- 
sumption of isotrop ic random motions, (e.g. JOnken et al.l2 004; 
IGreene & H0II2O O6) yielding: 



M BH = (9.7±0.5)xlO 



6 / 



^5100 



V 10 44 erg s" 



0.519±0.07 /t^wtt,,, \ 2 

/ FWHMh, 
^ 10 3 kms" 



(1) 

The masses calculated using this formalism are displayed in 
Table 3 and Figure la. Note that in the case of B 1422+231, Ha 
is outside of our observing window. In this case, we use the H/3 
line, and the following (very similar) formalism: 



M BH = (9.1±0.5)xl0 6 



^5100 



10 44 erg s" 



0.519±0.07 



/FWHM H/3 
^ 10 3 kms" 



(2) 

The H/3 masses are presented in Figure lb. 

The luminosities used here are derived from the HST pho- 
tometry of the CASTLES sample, as described in detail in P06. 
Briefly, the quasar and host galaxy images are optimized along 
with the lens m odel using LENSFIT (a variant of GALFIT; 
Peng et ai]|2002l) . Th e quas ars are modeled as point sources, 
the galaxies as Sersic ( 1968) functions, and the lensing masses 
as singular isothermal ellipsoids. The demagnified and de- 
blended V,I,H quasar magnitudes are modeled as a power-law, 
with the slopes presented in Table 2, and the luminosities come 



from this fit. We note that our spectroscopically derived slopes 
agree nicely with those from broad-band photometry. We find 
(a op t/apo6) = 1 ±0.5, where f\ oc A Q . The HST //-band images 
correspond closely to restframe V-band at z ~ 2, so the lumi- 
nosities are very insensitive to uncertainties in the powerlaw 
slope. 

We find a mean ratio of (log(M H a/Mpo 6 )) = -0.02 ±0.5 be- 
tween the two mass estimates. Although the scatter is large, 
we find no evidence for a systematic offset in the C IV-based 
masses. Therefore, to the extent that virial mass estimates 
have merit in this mass, luminosity, and redshift regime (as yet 
untested directly), our new spectroscopic observations confirm 
the results presented in P06. 

3.4. Uncertainties 

Uncertainties in line widths are difficult to estimate. In our 
case, the signal-to-noise (S/N) ratios of the spectra are not very 
high (see Table 1), which contributes substantially to the error 
budget. We use Monte Carlo simulations to estimate the mag- 
nitude of the uncertainties due to finite S/N. For each observed 
line, we generate 1000 mock spectra with the same shape as 
our best-fit model and the same noise as our observed spec- 
trum. In order to model residual uncertainties in continuum 
subtraction, we include a broadened Fe II spectrum with an am- 
plitude that scatters around 10% of the peak flux. Using the 
same procedures as in the real data, we then measure the best- 
fit parameters for each of these simulations. The uncertainty is 
determined from the distribution of simulated line dispersion to 
be half of the width encompassing 68% of the simulated galaxy 
measurements, and ranges from ~ 10-40% of the measured 
value. Now, in addition to noise, the treatment of narrow emis- 
sion lines can be a significant cause of systematic uncertainty. 
Thus we perform a second fit with the narrow component of 
each line turned off, and the difference between the two repre- 
sents a second estimate of the uncertainty. We adopt the larger 
of these as our final error in line width. 

At the same time, quasar variability adds additional un- 
certainty since the rest-frame optical and UV spectra ana- 
lyzed h ere were not obtained cont emporaneously. I Wilhite et al.l 
d2007l) and iDennev et ail d2009al) both find that variability in 
06±0 fme width tends to be small (~ 30%) and thus only con- 
tributes ^0.1 dex scatter to BH mass estimates. Another 
feature of quasar spectra that complicates emission-line-width 
measurements is the presence of blue-shifted absorption fea- 
tures with velocity widths of tens (narrow) to tens of thou- 
sands (broad) of km s . Both H1413+117 and (to a lesser 
extent) PG1 115+080 display broad absorption features, while 
Q0957+561 and HS08 10+2554 have significant narrow absorp- 
tion. In these cases we simply mask the regions from the fit, but 
note that particularly in the case of H1413+117, our ability to 
measure a reliable C IV width is seriously compromised, since 
even the line center is not well-defined. 

In addition to the line widths presented here, the BH masses 
depend on the luminosity of the quasar. We briefly review the 
arguments that the lens modeling does not add significant un- 
certainties to our results and direct the interested reader to Ap- 
pendix B of P06 for more details. The primary source of uncer- 
tainty in the models is that gravitational time delay and lensing 
substructures can lead to anomalous quasar magnification ra- 
tios, which translate into an uncertainty in the quasar luminos- 
ity. However, this error should be no larger than 0.1-0.2 mag, 
and thus does not contribute significantly to the BH mass un- 
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certainty. Actual measurement uncertainties from model fitting 
are tiny for the quasar. 

In Table 3 we present the formal uncertainties in BH mass 
arising from errors in the FWHM measurements and the for- 
mal uncertainty in the slope of the radius-luminosity relation. 
We wish to emphasize that the true errors in BH mass are 
probably dominated by systematic uncertainties arising from 
our ignorance of the structure and kinematics of the broad- 
line region that translate into errors in inferring both its true 
extent and velocity field based on the observations dKrolikl 
1200 It ICollin et al]l2006t iGreene & Holl2006l) . At present, we 
really have no concrete confirmation that (a) the same radius- 
luminosity relation applies to quasars at these luminosities (al- 
though see iKaspi et al.|[2Q07h nor that (b) it is meaningful to 
assume that the BLR gas is in virial equilibriu m and not, for 



instance, do minated by a massive outflow (e .g., Baldwin et al 
1996[lRichards et al.l2 002a: Leighl v & Moord2004tlFine et al 
20081) . 



4. DEMOGRAPHICS OF LENSED QUASARS 

To summarize, we have obtained Ha-based BH masses for a 
large fraction of the lensed quasars with measured host galaxy 
luminosities from P06. While the scatter between Ha-based 
and C IV- or Mg II-based BH masses is quite large (~ 0.5 dex), 
we do not see evidence for a systematic bias in mass. Thus, we 
do not alter the result of P06 that Mbh - ^buige relations appear 
to evolve with cosmic time. 

It is instructive to now compare the properties of the 
lensed sample with the general quasar population at 1 < 
z < 4. A large number of studies have looked at the 
distributions of B H mass and Eddington ratio of lumi- 
nous quasars (e.g. , iMcLure & Dunlod 120041: iKollmeier et ail 
200 61: iNetzer et al.1 120071 : iFine et alj 120081; IShen et alj 12008b; 



Gavignaud et al.1 120081) including some work on narrow - 
line quasars dGreene et al" | [20091 iLiu et al.1 12 009) as well 
as mass functions for quasars (e.g., Greene & Hoi [2007b; 



IVestergaard et al.ll2008t iKelly et al.ll2009l) . In most cases the 
BH masses are based on C IV and Mg II from observed optical 
spectra. Generally, all works find the same surprisingly nar- 
row range in measured line width, and correspondingly narrow 
distributions in M BH and iboi/^Edd at all redshifts (although, as 
Gavignaud et al. point out, the derived distribution is quite sen- 
sitive to the assumed slope in the radius-luminosity relation). 

For our purpose it is easiest to compare with the results of 
Netzer et al., as their study is also based on H(3 observations. 
For convenience we adopt their assumption that L5100 is ~ one 
seventh of the bolometric luminosity and we take the Edding- 
ton luminosity to be 1.26 x 10 38 (M BH /M Q ) erg s" 1 . Despite 
completely different selection criteria, the distributions in lu- 
minosity, and correspondingly inferred BH mass and Eddington 
ratio, are quite similar across the two samples (Fig. 8). For the 
Triplespec sample, the median (Mbh) ~ 2 x 10 9 M and the me- 
dian L5ioo~ 6 x 10 45 erg s" 1 yield a typical Eddington fraction 
of ~ 30%. If we take the entire P06 sample, using Ha-based 
masses when available, we find a median mass of ~ 10 9 Mq, a 
median L5100 of 2 x 10 45 erg s , and a median Eddington ratio 
of 10%. 

Unlike other studies, Netzer et al. actually argue that 
their sample displays a broad distribution in Eddington ra- 
tio when compared to theoretical expectations that luminous 
quasars should be most r eadily observed at or near their Ed- 
dington luminosities (e.g., Marc oni et~aLll2004t lMerlonil l2004; 



HopkinsetalJ[2006). We find the observed distributions un- 
expected for a somewhat different reason, namely the uni- 
formly high BH masses. Assuming that the Eddington limit 
strictly applies and given the lower flux limit of the lensed sam- 
ple we could detect BHs with masses as low as ~ 10 8 Mq. 
Furth ermore, our intuition from the local BH mass function 
(e.g., lYu & Tremaiiiell2002t iMarconi et ai]|2004l) suggests that 
10 8 Mq BHs ought to be far more numerous than 10 9 Mq sys- 
tems due to the exponential decline in the space density of the 
most luminous galaxies. Given that supercritical accretion may 
be observed locally (e.g.. [Pounds et alj|1995l: iMineshige et al.l 
l2000t iDesroches et al.l 120091) . we find the preponderance of 
~ 10 9 Mq BHs radiating at ~ 10% of their Eddington lumi- 
nosity to be noteworthy. 

We can be sli ghtly more quantitative. Following, e.g., 
ISomerviiiel d2009l) . we can transform th e observed z ~ 2 galaxy 
mass function (e.g.. iFontana et aLll2006l) into the expected BH 
mass function. For reference, at the median luminosity of 
the Triplespec sample, the Eddington limit sets a bound of 
~ 4 x 10 8 Mq on observable BHs. If the relation between 
Mbh and galaxy ma ss were identical to what it is today (e.g., 
Gtilt ekin et al.ll2009l) . then the BHs with mass w 4 x 10 8 M Q 
ought to be ~ 50 times more common than the median ob- 
served mass of M B h~ 2 x 10 9 Mq. Note that we are con- 
servatively adopting the median rather than the minimum ob- 
served luminosity and assuming that the Eddington luminosity 
applies strictly. Taking the observed masses at face value for 
the moment, and assuming that the selection of lensed quasars 
is simply a random selection of optically luminous quasars at 
1 < z < 4, then we infer that either BHs radiating at 30% of Ed- 
dington are ~ 50 times more numerous than Eddington-limited 
objects or that the zeropoint in the Mbh _ ^buige relation has 
evolved by a factor of ~ 3 to the present day. This latter, of 
course, is the suggestion made by P06. 

Obviously, such constraints are not particularly stringent at 
the moment. For one thing, the selection of lensed quasars is 
complicated to model, involving radio selection in some cases 
(whi ch will tend to bias samples toward more m assive systems; 
e.g.. iHeck man 1983; Mandelbau m et al.ll2009l) . Furthermore, 
as pointed out by Somerville, joint constraints from galaxy and 
quasar luminosity functions on their own cannot distinguish 
between zeropoint evolution or increased scatter in BH-bulge 
relations at high redshift. Finally, we do not know that the 
quasar duty cycle is independent of mass. Nevertheless, we 
cannot help but wonder whether the narrow range in observed 
line width, and the unexpectedly high average BH mass, do not 
instead indicate a problem in virial mass estimators at high lu- 
minosity that may be resolved with better understanding of the 
physics of broad-line quasars. 

Ultimately, the question is whether or not the broad emis- 
sion lines in luminous quasars are dominated by virial motions. 
Of course, there is evidence for a non-virialized component i n 
the C IV line (e.g jBaldwin et alJll996l: iRichards et al.ll2002al) . 
On the other hand, there is clear evidence for virialized mo- 
tion in the broad-line regions of at least a few nearby lower- 
luminosity active galaxies (e.g., Peterson & Wandell [l999; 
lOnken & Peterson! 120021: iPeterson et al.ll2004l) . Als o, the local 
virial m asses seem to c orrelate with both a* (e.g., Shen et al.l 
2008a|) and L bu i g e (e.g., iKim et alj|2008bt iBentz et al.l l2009b). 
Unfortunately, similar analysis is not yet available in the lu- 
minosity range of interest to us. At present, all we can say is 
that the C IV-based masses alone are not causing a net bias in 
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Fig. 7. — (a): We compare BH mass estimates based on the Ha transition with those from P06. The P06 masses are based on C iv with one exception (red 
open circle; Mg n) while our masses are based on Ho with one exception (magenta open square; H/3). The unity relation (solid line) is shown to guide the eye. 
While the two mass estimates are not strongly correlated, we do not find any evidence for a systematic offset between the two (the median Mn a /Mpo^ = 1.0 ± 0.4). 
H1413+1 17, due to its broad-absorption system, is highlighted as a filled symbol, (b): Same as (a), except using H/3 as the virial indicator. In this case objects with 
only partial observations of the H/3 line are highlighted with crosses. As above the P06 mass that is based on Mg u is identified with a red circle. We find a median 
M Hl3 /M Pa6 = 1.0±0.5. The filled triangle is H1413+117. 
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FIG. 8. — (a): We reproduce Figure 2 from Netzer et al. (2007, black dots; based on H/3) and include the targets from this study (large black open squares), as 
well as the entire P06 sample, differentiating between Mg n-based (small red circles) and C iv-based (small blue triangles) masses. Here we plot monochromatic 
5100 A luminosity versus the Eddington ratio. We follow Netzer et al. and assume that the bolometric luminosity is seven times £5100- Both our mass estimates and 
those from Netzer et al. are based on Balmer emission lines (H/3 in their case), (b): Mbh versus Lbol/^-Edci; symbols as in (a). We see that the distributions in each 
diagram are similar between the two samples, although the lensed systems extend to objects that are factors of two to three times fainter. 
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the virial masses relative to the Balmer lines. True evolution 
in BH-bulge scaling relations is by no means certain; far more 
pernicious sources of uncertainty remain, including potential 
biases in sample selection. At the minimum, to mitigate these 
concerns, quasar samples with identical selection at multiple 
redshifts are needed. 

In closing, we note that BH-bulge relations are not the only 
ones purported to display unexpected evolution since a red- 
shift z»2. Recent work has shown that, at a fixed mass, el- 
liptical galaxies were a factor of two to five smaller at red- 
shifts of one and t wo respectively than they are today (e.g. , 
Truiillo et alj|2006tlvan Dokkum et al.ll2008l:lFranx et alJ 2008 ; 
van der Wei et al.l2 008; Da mianov et alj2009l) . At first it seems 



counterintuitive that the most massive elliptical galaxies, which 
seem to have fo rmed the bulk of the ir stars rapidly at an ear- 
lier epoch (e.g., Thomas et al] |2005l) . should grow less dense 
with time. However, a large number of minor mergers can 
very efficiently build the outskirts of elliptical ga l axies (e.g., 
Bovl an-Kolchin & Mall2007t iBezanson etai]|2009t iNaab et al.l 



2009). Interestingly, these gas-free minor mergers would need 
to grow the total galaxy mass by fact ors of two to three in 
order to match local observations (e.g.. IBezanson et al1 l2009). 
This is possible under a scenario where the BH-to-TOTAL stel- 
lar mass relation is steeper than linearit y in low mass galax- 
ies, which is obser ved in nearby galaxies (Giilt ekin et al.l 2009; 
iGreene et al.l2008l) . In gas-free merging there is no correspond- 
ing BH growth via accretion, which coul d facilitate a boost i n 
the BH-bulge ratio at late times (see also [Hopkins et al. 2009). 
As argued bv lPengl (12001 . numerous minor mergers will tend 
to drive the ratio of BH to galaxy mass towards a value of unity 
slope (e .g., his Fig. 2 b a nd 4a). Alternativel y, as suggested 
by bofh lCrotonl d2006l) and iJahnke et al.1 d2009l) . it may be that 
disk components are common in the high redshift galaxies but 
are subsequently subsumed into the main elliptical galaxy body. 
We mention this apparent coincidence in passing because it is 
intriguing, although we note that a variety of outstanding uncer- 
tainties are yet to be explored as far as the structural evolution 



of elliptical galaxies is concerned. 

5. SUMMARY 

We revisit the BH mass estimates for a sample of 
lensed qua s ars with high-fidelity host-galaxy luminosities from 
Peng_et_alJ ([2006b). While the Balmer-based masses presented 
here are arguably more robust than the UV-based estimates, we 
find no evidence for a systematic difference in BH masses based 
on the two methods. If we can take the Balmer- based virial 
masses at face value, then we confirm the result of iPeng et al.l 
(2006b) that BHs appear to be overly massive relative to their 
hosts at high redshift. Intriguingly, the minor mergers that 
are invoked to puff up compact elliptical galaxies at late times 
would be very effective at boosting the galaxy to BH mass ratio 
as well. However, as pointed out many times, the persistently 
narrow range in observed line width for luminous quasars at 
high redshift provides substantial cause for concern in the ve- 
racity of the virial masses. Furthermore, it is not clear that 
we have yet assembled consistent comparison quasar samples 
across cosmic time with which to compare the BH-bulge scal- 
ing relations. Our work is a necessary, but not sufficient, step 
in determining the true evolution of BH-bulge relations with 
cosmic time. 

The referee provided a very careful and insightful reading 
that substantively improved the manuscript. Both J. E. G. and 
C. Y. P. are grateful for many stimulating conversations with 
L. C. Ho and we thank A. J. Barth and C. A. Onken for very 
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sioning data for this project and for invaluable help with the re- 
duction software. We thank H. Netzer for kindly providing his 
data table for Figure 8. Support for J. E. G. was partially pro- 
vided by NASA through Hubble Fellowship grant HF-01196 
awarded by the Space Telescope Science Institute, which is op- 
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APPENDIX 



Triplespec Spectra 

In Figure Al(a-d), we present the rest of the Triplespec spectra for completeness. Note that the Triplespec spectra for SBS0909 
and HS0810 are shown in Figure 1. 
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FIG. 9. — Fits to the continuum (left), H/3 (middle), and Ha (right) lines from the Triplespec data that do not appear in Figure 1. We show the data (solid 
histogram), the total model (thin solid), the broad- (dashed) and narrow-line (dotted) model components, and residuals below (thin solid histogram). Data are 
plotted with an arbitrary scale in fx- Spectral regions that are masked in the fit do not appear in the residuals (e.g., the red wing of H/3 for Q0142-100). Also note 
that in a couple of cases the [O m] fit is based exclusively on the A4959 line; no redshifts are derived for these targets. 



Table 3. Line- width Measurements 



Name 


Z FWHM[oni] 


FWHMho 


FWHM H/ 3 


FWHM Mg ii 


FWHMcni] 


FWHMciv 


^£5100 


log M B H.P06 


log Mbh.Hc log t 


(1) 


(2) (3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


(11) ( 



PMNJ0134-0931 


2.22 


1410 


5.1±0.5 


4.2± 1.4 
















Q0142-100 


2.7 




3.8± 0.3 


2.7± 0.6 




6.0± 0.2 


6.5± 1.3 


46.42 


9.4 


9.4± 0.3 


9.1 


SDSS0246-0825 


1.69 


540 


2.5± 0.2 


2.5± 0.6 


2.9± 0.1 


4.5±2.8 


5.0± 0.2 










MG0414+0534 


2.6 




5.3± 0.3 


10.8± 2.3 








45.73 


9.3 


9.4± 0.3 




HS0810+2554 


1.51 


490 


3.8± 0.0 


4.4± 1.1 


3.6±0.1 


4.6± 0.4 


3.9± 0.1 










SBS0909+523 


1.38 


410 


3.1±0.6 


4.1± 1.6 


3.5±0.5 


4.1± 1.1 




45.96 


9.6 


9.0± 0.3 


9.2 


FBQ095 1+2635 


1.25 


290 


6.3± 0.4 


7.6± 0.8 


6.0± 0.3 


7.2± 0.8 




45.48 


8.9 


9.4± 0.3 


9.5 


Q0957+561 


1.41 


570 


3.0± 0.2 


3.3±0.9 


3.3± 0.2 


5.2± 0.5 




46.25 


9.3 


9.1±0.3 


9.2 


J1004+1229 


2.7 




3.8±0.1 


5.9± 1.5 








46.12 


9.3 


9.3± 0.3 


9.6 


HE1 104-1805 


2.32 


970 


4.7± 0.2 


3.8± 0.9 








45.59 


9.4 


9.2± 0.3 


8.9 


PG1 115+080 


1.7 




4.0± 0.1 


4.4± 0.2 


4.0± 0.2 


5.7± 3.1 


6.7± 1.0 


45.15 


9.0 


8.8± 0.2 


8.8 


B 1152+200 


1.0 




7.1±2.6 








3.3±0.6 










H1413+117 


2.55 


1650 


5.3± 0.8 


6.7± 1.9 




8.0±3.2 


3.3± 0.6 


45.67 


8.4 


9.3± 0.3 


9.5 


B 1422+231 


3.63 


2720 




6.1±2.2 








46.65 


9.7 




9.9 


FBQ1633+3134 


1.52 


1720 


4.1±0.7 


4.6± 0.9 


4.1 ±0.5 


6.4± 0.6 


3.8± 0.7 


45.81 


9.2 


9.2± 0.3 


9.2 


Q2237+030 


1.70 


1070 


4.8± 0.6 


3.8± 1.4 

















Note. — Col. (1): Name. Col. (2): Redshift. Cases with three significant digits were measured from [O III] while cases with two come from the Castles 
(http://www.cfa.harvard.edu/castles/). See P06 for references. Col. (3): FWHM of [O III] A5007 A line (km s _1 ), with the instrumental resolution of <r» £b 36 km 
traded in quadrature. Col. (4): FWHMn a (10 3 km s -1 ) measured from the Triplespec data. Col. (5): FWHMh/3 (10 3 km s -1 ) measured from the Triplespec data. Asterisi 
objects with incomplete spectral coverage of the H/3 line. Col. (6): FWHMMgii (10 3 km s -1 ) measured from the SDSS spectrum. Col. (7): FWHM C ni] (10 3 km s -1 ) i 
from the SDSS spectrum. Col. (8): FWHMciv (10 3 km s -1 ) measured from the SDSS spectrum. Col. (9): Monochromatic luminosity (erg s _1 ) at A = 5100 A as meast 
photometry (P06). Col. (10): Log {M^IMq) from P06; typically based on the C IV line. We do not quote errors in the BH masses, since systematic uncertainties domi 
typically quote a factor of four error as derived by Vestergaard & Peterson 2006. Col. (11): Log (Mbu/Mq) calculated using Equation 1 and the FWHMhc* from Col. (4) c 
with the luminosity from Col. (9). Note that the errors in the Ha masses shown here include only measurement uncertainties. Col. (12): Log (M^IMq) based on H/3. In 
Equation 2 and FWHMh^ from Col. (5) are combined with the luminosity from Col. (9). 



